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Abstract 

In this paper we analyze the dynamic stability of the hovering 
magnetic top from first principles without using any preliminary as- 
sumptions. We write down the equations of motion for all six degrees 
of freedom and solve them analytically around the equilibrium solu- 
tion. Using this solution we then find conditions which the height of 
the hovering top above the base, its total mass, and its spinning speed 
have to satisfy for stable hovering. 

The calculation presented in this paper can be used as a guide to 
the analysis and synthesis of magnetic traps for neutral particles. 



*Also with the Department of Physics, University of California, San Diego, La Jolla, 
92093 CA, USA. 
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1 Introduction 



1.1 The hovering magnetic top. 

The hovering magnetic top is an ingenious device that hovers in mid-air while 
spinning. It is marketed as a kit in the U.S.A. and Europe under the trade 
name Levitron^^- [j], |2j and in Japan under the trade name U-CAS0. 
The whole kit consists of three main parts: A magnetized top which weighs 
about 18gr, a thin (lifting) plastic plate and a magnetized square base plate 
(base). To operate the top one should set it spinning on the plastic plate 
that covers the base. The plastic plate is then raised slowly with the top 
until a point is reached in which the top leaves the plate and spins in mid- 
air above the base for about 2 minutes. The hovering height of the top is 
approximately 3cm above the surface of the base whose dimensions are about 
lOcmx 10cmx2cm. The kit comes with extra brass and plastic fine tuning 
weights, as the apparatus is very sensitive to the weight of the top. It also 
comes with two wedges to balance the base horizontally. 



1.2 Qualitative Description. 

The physical principle underlying the operation of the hovering magnetic top 
relies on the so-called 'adiabatic approximation' f|, ||, |], [7[: As the top is 
launched, its magnetization points antiparallel to the magnetization of the 
base in order to supply the repulsive magnetic force which will act against the 
gravitational pull. As the top hovers, it experiences lateral oscillations which 
are slow ( ^i a ^ era i — 1Hz) compared to its precession (^p recess i on ~ 5Hz). 
The latter, itself, is small compared to the top's spin (^ S pj n ~ 25Hz). Since 
^spin ^ ^precession ^he top * s considered 'fast' and acts like a classical 
spin. Furthermore, as ^p recess i n ^ ^lateral ^ ms s P m ma Y De considered 
as experiencing a slowly rotating magnetic field. Under these circumstances 
the spin precesses around the local direction of the field (adiabatic approx- 
imation) and, on the average, its magnetization \x points antiparallel to the 
local magnetic field lines. In view of this discussion, the magnetic interac- 
tion energy which is normally given by —y, ■ H is now given approximately 
by \i \H\. Thus, the overall effective energy seen by the top is 

E e ff ~ mgz + fi \H\ . (1) 
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By virtue of the adiabatic approximation, two of the three rotational degrees 
of freedom are coupled to the transverse translational degrees of freedom, 
and as a result the rotation of the axis of the top is already incorporated in 
Eq.(JJ). Thus, under the adiabatic approximation, the top may be considered 
as a point-like particle whose only degrees of freedom are translational. The 
important point of this discussion is the following: The energy expression 
written above possesses a minimum for certain values of fi/m. Thus, when 
the mass is properly tuned, the apparatus acts as a trap, and stable hovering 
becomes possible. 

As mentioned above, the adiabatic approximation holds whenever fi S pj n 3> 

^precession ana - ^precession ^lateral' ^ s ^precession * s mverse ly pro- 
portional to ^ S pj n , these two inequalities can be satisfied simultaneously pro- 
vided that the top is spun fast enough, to get ^ S pj n 3> ^precession; Du ^ no ^ 
too fast, for then ^precession > ^lateral" ^ ne reason f° r the lower bound on 
the spin is obvious: If the top is spun too slowly, then fi S pi n < ^precession 
and the top becomes unstable against rotations. The top will then flip over 
and will be pulled quickly to the base. This is the instability that is well 
known from classical top physics[[iT], I2| . The reason for an upper bound 



on the spin is quite different: If the top is spun too fast, the axis of the 
top becomes too rigid and cannot respond fast enough to the changes of the 
direction of the magnetic field. It is then considered as fixed in space and, 
according to Earnshaw's theorem|TU[ , becomes unstable against translations. 



1.3 The purpose and structure of this paper. 

The major problem with the adiabatic approximation is that it cannot predict 
exactly the allowed range of ^gpj n for which stable hovering occurs. It simply 
gives us estimates for this range, and for design purposes this may not be 
enough. The purpose of this paper is to give a quantitative description 
of the physics of the hovering magnetic top while it is in mid-air without 
using any preliminary assumptions (such as the adiabatic approximation). 
To do this, we first expand the magnetic field around the equilibrium point 
to second order in the spatial coordinates. Using the gradients of this field 
we then find the force and torque on the top and write down the equations of 
motion for all 6 degrees of freedom in vectorial form|8|]. Next, we solve these 
equations for the stationary solution and for a small perturbation around 
this stationary solution and arrive at a secular equation for the frequencies 
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of the various possible modes. The possible eigenmodes are either oscillatory 
(corresponding to stable solution) or exponential (unstable solution). The 
secular equation comes out linear in fi S pj n , and it is not difficult to write 
analytic expressions for the upper and lower bounds of fi S pj n - 

The structure of this paper is as follows: In Sec.@ we first describe our 
model and define our notations. Next, we derive the equations of motion 
for the translational and rotational degrees of freedom, and finally we solve 
these equations around the equilibrium position. In Sec.| we apply our results 
to the case of a disk-like top of radius a = 0.25R hovering above a circular 
current loop of radius R. We fix the equilibrium hovering height and plot the 
various frequencies of the stable modes as a function of ^ S pi n - We identify the 
origin of the various modes and comment on the way these modes are coupled 
to produce the minimum and maximum speeds at which the system becomes 
unstable. In particular the connection of these results with the adiabatic 
approximation will be discussed. Then we change the equilibrium height 
and note how the allowed range of ^ SD i n is affected. We shall show that 
above a particular height, where according to the adiabatic approximation 
hovering is not possible, hovering is still possible in a certain range, and 
that, furthermore, the lower bound on fi S pj n originates from a new mode 
coupling, not predicted by the adiabatic approximation. In the last section 
we summarize our results and discuss the possible uses of the derivation 
presented in this paper to the study of magnetic traps for neutral particle. 

2 Derivation and Solution of the Equations 
of Motion. 

2.1 Description of the model and notation. 

In this paper we analyze the dynamics of the top while in mid-air as is shown 
in Fig.^. We consider a symmetric top with mass m and a magnetic moment 
/i, rotating about its principal axis with angular speed Q s . The magnetic 
moment of the top is assumed to be concentrated at the center of gravity 
of the top and to point along the symmetry axis of the top. The latter is 
denoted by the unit vector n. The moment of inertia of the top about n is 
^3 whereas the secondary moment of inertia is I\. 

The spatial position of the top is given with respect to its equilibrium 
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position (point Q in the figure). Thus, z is the vertical displacement of the 
top and p is its radial displacement. We denote by z and p the unit vectors 
in the vertical and radial direction, respectively. 

Though this is not mandatory for our calculations, we assume for sim- 
plicity, that the magnetic field possesses cylindrical symmetry around the 
vertical axis. In addition, the direction of gravity is assumed to coincide 
with the symmetry axis of the magnetic field and to point downward. We 
denote by g the free fall acceleration. 



2.2 The equations of motion. 

As the symmetry axis of the magnetic field coincides with the direction of 
gravity, equilibrium is possible only along this axis. We therefore express the 
magnetic field, to second order in p and z, at the vicinity of the equilibrium 
position in terms of its derivatives along the z direction at the equilibrium 
position. This is possible due to cylindrical symmetry and to the fact that 
the Cartesian components of the field are harmonics. The result is 



H(p, z)= -\p (H' + H"z) p + [H + H'z + \H" (z 2 - ip 2 ) 



z. (2) 



Here H , H' and H" are the vertical magnetic field, its first and second 
derivatives along the z direction, respectively, at the equilibrium position, 
i.e. at point Q. 

The potential energy of the top is the sum of the magnetic interaction 
energy of a dipole with a field plus a gravitational term, i.e. 

E = —ph ■ H (p, z) + mgz. 

Consequently the force on the top is 

F = -VE = pS7 (h ■ H) - mgz, (3) 

whereas the torque is 

T = pnx H. (4) 

The equation for r =pp+zz, the radius vector of the center of mass of the 
top with respect to point Q, is given by 

d 2 r 

m-j^=pV(n ■ El)— mgz. (5) 
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To write a vectorial equation of motion (see for example ||) for the an- 
gular momentum, L , we note that it has two components in perpendicular 
directions. The first component is due to the rotation of the top around 
the n direction and is given simply by L n = I^QgU. The second component 
of the angular momentum is contributed by the change in the direction of 
the principal axis from n to n+dh. Since, by definition, n is a unit vector, 
it must point perpendicular to dn. Thus, Q± = \dn/dt\ is the angular ve- 
locity associated with the change of n. Since the direction of flj_ must be 
perpendicular to both dn and n we form the cross product fl± = hx (dn/dt) 
which incorporates both the correct value and the right direction. Multiply- 
ing Q± by I\ yields the second component of the orbital angular momentum, 
Lj. = hnx (dh/dt). Thus, L = L n + L± = I^VL s h+Iih x (dh/dt). Using this 
expression together with Eq.(P for the torque we find that the equation of 
motion for the angular momentum is 



dL _ d 
dt dt 



T „ A r „ dn 
hilsn+hnx — 
dt 



fih x H. (6) 



Eqs.|5],|6| together form a coupled system of equations for all 6 degrees of 
freedom of the top. In these equations r, n and Q s are the dynamical variables 
(a total of 6 degrees of freedom), H is the magnetic field which itself is a 
function of r, whereas m, /x, I\ and J3 are external parameters. 



2.3 Solution for the equations of motion 

The stationary solution of the problem is obvious: When the top is on the 
symmetry axis with its principal axis parallel to the z axis, two forces act 
on the top. The first force is the downward gravitational force mg and the 
second one is the upward magnetic force supplied by the external magnetic 
field, i.e., [iH'. Since these forces are colinear, no torque is exerted. We 
therefore look for a solution of the form 

p(t) = z(t) = 

n(t) = z (7) 
Q s (t) = Q s = Const. 

Inserting this solution into Eqs.|,^| yields identities for 5 degrees of free- 
dom. The equation for z, on the other hand, gives the expected equilibrium 
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condition 

pH' = mg. (8) 

To investigate the stability of the stationary solution, we now add first 
order perturbations to the equations of motion. Thus, we make the following 
substitutions 

z(t) -> + 5z{t) 
p{t) - + 5p(t) 

n s (t) -+ n s (t) + sn(t) { ! 

where, since n is a unit vector, 5n must be perpendicular to z. Substituting 
Eq.[| into Eqs.|5],|6| and expanding to first order in the perturbations gives 

d 2 5z 

= ^ H "5z (11) 
dt 2 ^ v ; 

£-0 

l3 ^ s ~dT + X ~d# = ~^ H ^ x Sfl ~~ 2^'^ X ^ 

We note that (at least to lowest order) the motions in the z direction and 
the rotation around the n axis are decoupled from the other degrees of free- 
dom. Furthermore, according to Eq.[ll| the motion in the Indirection is stable 
provided that H" < 0. We now focus our attention on the remaining four de- 
grees of freedom: Note that the right-hand side of Eq.[l^ contains two terms: 
The first is the adiabatic term which tends to stabilize the top against lateral 
translations by tilting the axis of the top. The second term, which we call 
Earnshaw's term, tends to destabilize the top and to take it away from the 
equilibrium position. Solving Eq.|10| for 5h and substituting it into Eq.|13| 
results in a fourth order equation for the radius 5p\ 

2ml 1 z x <5p< 4 ) + 2mI 3 Q s 5^ + 
+p [hH" + 2mH] z x 5jP^ + pI 3 n s H"5^ . (14) 
-p 2 \l (H') 2 - HH"} zx5o = 
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The possible solutions of this equation are linear combinations of a steady 
rotation of Sp around z at angular velocity u (this is possible because of the 
cylindrical symmetry of the field; otherwise, one should write two equations 
for the two components of Sp). We thus set 



dSp 



uoz x Sp. 



(15) 



Substitution of Eq.[L5| into Eq.[14] results in the following secular equation for 
the eigenfrequencies: 



-p [hH" + 2mH] lu 2 + pI 3 tt s H"uj 



\ (H'f - HH" 



(16) 



0. 



This fourth-order equation for u has four real roots (or eigenfrequencies) 
whenever the system is stable. Since we are looking for the range of Q s for 
which the system is stable, we take another point of view and express fl s in 
terms of u. The resulting equation after using mg = pH' is 



2hH'uj A - g [hH" + 2mH] a; 2 



mg 



hH' 



HH' : 



H' 



2I 3 H'lu 3 - gI 3 H"uj 



(17) 



3 Application to a disk-like top above a cir- 
cular current loop. 

As an example we now apply Eq.[l7] to the case of a disk-like top of radius 
a. Consequently, /3 = ma 2 /2 and I\ = ma 2 / 4. The source of the magnetic 
field is taken as a horizontal current loop (or alternatively, a vertically uni- 
formly magnetized thin disk) of radius R. The vertical magnetic field and its 
derivatives along the axis at a height h above the loop are therefore given by 



H = -H l + (h/R) 



-3/2 



H' 



Mo 

R 



H" 



Mo 

R 2 



(h/R) 
4 (h/R) 



(h/Rf 
"l + (h/Rf 



-5/2 



(18) 



-7/2 
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Note that the magnetic field was chosen to point downward. Thus, in order to 
get an upward repulsive magnetic force, the magnetic moment should point 
upward. The sign convention is chosen such that h and \i are positive, H is 
negative and, consequently, H' is positive. For stability in the z direction we 
require H" < 0. This occurs, according to Eq.|l8|, as long as h > h min = 0.5R, 
and sets a lower bound for the height at which stable hovering is possible. 
Taking h = 0.55R and a = 0.25R (these are approximately the parameters 
for the Levitron) inside Eq.|17] and plotting Q s versus u> yields the solid line 
shown in Fig.^. Note that both f2 s and 10 are normalized to ujq = yJg/R. 
From this figure we learn that (for h = 0.55R and a = 0.25R) whenever 

7.lJg~/R<n s < 19.3 Jg~/R 

there are four real solutions u which correspond to four stable modes. The 
frequency of the fastest mode goes asymptotically to the dotted line Q s = 
(Ii/I 3 )uj as Q s ^> uq. This is nothing but the 'fast precession' rotational 
mode encountered in classical top physics [[iT| . The frequency of the next 
fastest mode goes roughly like u ~ fiH/I 3 Q s (compare it to the dash-dash- 
dotted line). This is the well known 'slow precession' rotational mode of a 
classical topjTTJ. The coupling between these two modes produces the min- 
imum speed for stability £\ m i n = 1 X^jgjR. For comparison purposes we 
have also included in our plot (dashed line) the resultant mode frequencies 
when we set H' and H" = inside Eq.([l6|). This is equivalent to solving the 
problem of a magnetized top, in a homogeneous magnetic field. In fact, these 
mode frequencies may be obtained as the roots of Q s = (h/h) uj + fiH/I 3 u. 
The minimum value of this expression occurs for u 2 = fiH/Ii and is given 

by fi s ,min = \J ^HIi/I 2 , which is the minimum speed of a classical top in a 
homogeneous field. By comparing the mode frequencies for the homogeneous 
field to the mode frequencies for the inhomogeneous field we see that (as far 
as the minimum speed is concerned) the minimum speed in both graphs is 
almost the same. The two slowest modes in the solid-line plot are the two 
vibrational modes of the top. It is clearly seen that one of them is strongly 
coupled to the slow precession mode. This coupling is responsible for pro- 
ducing the maximum speed for stability Q s ,max = lQ.Sjg/R, as was already 
predicted by the adiabatic approximation. But unlike in the adiabatic ap- 
proximation, we now have a way to find analytically both the minimum and 
maximum speed in a single stroke. 
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Our next step is to study how the allowed range of f2 s depends on the 
equilibrium height h. In Fig.[| we have plotted the mode frequencies versus 
n s for h = 0.55R, O.QR, 0.Q325R , 0.64R and 0.65.R. In Fig.|we have plotted 
the minimum and maximum speeds for each equilibrium height starting at 
h = 0.5R. Recall that stability along the z direction requires that h > 0.5R. 
Thus, the two curves plotted in Fig.|] together with the h = 0.5R line define 
a closed region in the Q s -h plane. Each point inside this region corresponds 
to a stable hovering solution whereas each point outside this region belongs 
to an unstable solution. From this figure we may deduce that the range of 
heights for which stable hovering is possible is given by 



The figure also shows that as h increases above its minimum value ^ m j n = 
0.5R, the minimum allowed speed ^ sm j n i decreases slightly, whereas the 
maximum allowed speed fl s ,max increases. As h increases further, the cou- 
pling between the two vibrational modes becomes stronger, and when h ex- 



ceeds a first critical value h Ci \ = y2/5R, this coupling gives rise to the ap- 
pearance of a new minimum ^ sm j n 2 of £l s (u), which increases steeply with 
increasing h. As the height exceeds the slightly larger value h c> 2 = 0.634i?, 
this new minimum Q s m j n 2 becomes higher than the minimum speed fl s m [ n \ 
determined by the coupling between the rotational modes, and thus limits 
the stability of the top. 

At h — h Cj 3 = 0.658.R, the minimum speed ^ sm i n 2(^) crosses the maxi- 
mum speed fi s ,max(h), such that stable hovering is not possible for h > h Ct3 
and Q s > ^mir^lX^)- 

It is important to note that the adiabatic approximation does not predict 
this new coupling. According to the adiabatic approximation the system will 
be stable against lateral translations whenever the curvature of the effective 
energy (given in Eq{l]) along the p direction is positive. It can be shown that 



this is satisfied whenever h < J2/5R. Thus, the adiabatic approximation 



predicts that hovering above h = y2/5R is not possible. In the present 
calculation we find, however, that the top is also stable above this height. 
This is due to the splitting of the vibration degenracy by the coupling to 
the precessional mode present when the spin is finite. This has already been 
pointed out by M. V. Berry § who treated it in terms the phenomenon 
of geometrical magnetism. Since each equilibrium height corresponds to a 



0.5R < h < 0.658R 
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different repulsive magnetic force (and hence to a different mass) we can, 
alternatively, specify the allowed range of mass required for stable hovering 
by using Eq.|8| and Eq.|IS] (note that the maximum height corresponds to the 
lightest mass and the minimum height corresponds to the heaviest mass). 
The results are plotted in Fig. [|. From the figure we learn that, for a = 
0.25.R, the mass of the top must be such that 

0.803m < m < 0.859m 

where 

fiH 

This corresponds to a mass tolerance of 

Am mmax ~ ^min ^ 

= 100— Sf£ mm _ 7% ig 

m m max 

Experimentally the tolerance is only about 1% which seems like a large dis- 
crepancy. Note, however, that the lower mass region (of less than 0.818mo) is 
difficult to access. Also, the top retains a finite kinetic energy when launched 



into the trap which, as has been discussed by us^3|, further decreases the 
mass tolerance in practice. Last but not least, the theoretical mass tolerance, 
given in Eq.([l9|) decreases drastically with the tilt of the base, and goes to 
zero for a tilt of about 0.45 degrees. This, in turn, makes it difficult to real- 
ize, in practice, the above theoretical mass tolerance, even with the leveling 
wedges supplied with the Levitron. 

Note also that the temperature coefficient of the magnetization of the 
ceramic magnet from which the Levitron is made is about 0.2% per degree 
Celsius. As the magnetic force varies with the square of the magnetization, 
temperature changes may easily 'throw' the mass out of range. 

To overcome these problems the kit comes with a set of light washers to 
tune the mass properly. 



4 Summary. 

In this paper we have analyzed the hovering magnetic top while it is in mid- 
air without using any preliminary assumptions. To do this we expanded the 
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magnetic field around the equilibrium point to second order in the spatial 
coordinates. Using the gradients of this field we then found the force and 
torque on the top and wrote down vectorial equations of motion for all 6 
degrees of freedom. Next, we solved these equations analytically for the sta- 
tionary solution and for a small perturbation around this stationary solution 
and arrived at a secular equation for the frequencies of the various possible 
modes. We then applied the solution to the case of a disk-like top hovering 
above a circular current loop, and were able to predict both the minimum 
and maximum allowed speeds for stable hovering to occur. 

Although the numerical results we have presented in this paper refer to 
the case of a disk-like top (/3//1 = 2), our theory is valid for the whole 
range of the anisotropy parameter J3//1 of the symmetric top, and is there- 
fore more general than the model studied in Ref.|| who approximated the 
top by a classical spin. We recover their results by setting I\ = in our 



calculations [13]. Also, our analysis yields naturally the minimum spin for 
stability, which is zero for the I\ = case, and thus is put in 'by hand' 
in the previous || treatment. The determination of the maximum spin, as 
a function of the anisotropy parameter I3/I1, is of particular interest for 
cigar-like tops {I3/I1 < 1), as for such tops the minimum speed behaves 
drastically different than for disk-like tops (2 > I3/I1 > 1), as we have shown 



elsewhere HR FU; 



An interesting corollary of our analysis is a limitation on the radius of 
the disk-like top. One would have expected that increasing the radius of 
the top will facilitate the operation of the Levitron as the moment of inertia 
increases, thus reducing the minimum speed. However, the maximum speed 
decreases faster when the radius increases, with the result that 

^s,max ,, „„ R 
~ 0.66—. 

^s,min r 

Thus, the top cannot exceed two thirds of the base. Another detrimental 
effect of increasing the radius is the reduction of the hovering time||17|| due 



to the increased effect of air friction. Moreover, air friction also modifies 
the stability analysis: So far we carried out a detailed study of the effect of 
friction on stability only for the I\ = model. In this case we found[|16| that 



with friction the top is always unstable even if the friction is infinitesimaly 
small and even if we invoke translational viscosity only. We expect this 
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behavior to occur also for /i ^ indicating that spin traps are different in 
character from potential traps in which friction 'increases' stability. 

We have disregarded here the angular momentum carried by the electrons 
responsible for the ferromagnetic moment of the top. Although this is very 
small compared to the orbital angular momentum of the Levitron, the two 
angular momenta may become comparable for sufficiently small tops, which 
would result in a left-right asymmetry. Furthermore, the possibility of levi- 
tation based only on the electronic angular momentum arises. This will be 
discussed elsewhere fi~6[ . 

Also, our treatment is completely classical. As size decreases, quantum- 
mechanical effects may become important. In particular, a sufficiently small 
particle, for example an atom, will only have a finite life-time in the trap due 
to quantum-mechanical effects, which will be considered elsewhere [[161. 
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Figure 1: The hovering magnetic top near its equilibrium position (point Q) 
and notations. 



15 




Figure 2: Mode frequencies vs. spin for a disk-like top of radius a hovering 
at a height h above a circular current loop of radius R. Solid line: inhomo- 
geneous field, dashed line: homogeneous field, dotted line: fast precession 
limit, dashed-dotted line: slow precession limit. 
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Figure 3: Mode frequencies vs. spin for a disk-like top of radius a hovering 
at a height /i = 0.55/?, 0.6R, 0.6325i?, 0.64/? and 0.65.R above a circular 
current loop of radius R. 
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Figure 4: Boundaries of the stability region in the Q s -h plane. 



18 



0.88 
0.87 
0.86 
0.85 
0.84 
0.83 
0.82 
0.81 
0.80 









m = 0.859m 

max 



maximum mass 

• maximum spin 

• minimum spin 



m a = \mjgR 
a = 0.25R 

®=g/R 




10 12 14 16 18 20 22 



£1/0) 



Figure 5: Boundaries of the stability region in the Q s -h plane. 
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